hep-ph/0603077 

Little Higgs Dark Matter 



Andreas Birkedal^, Andrew Noble^, Maxim Perelstein^, 
and Andrew Spray^ 

^ SCIPP, University of California, Santa Cruz, CA 95064 
Cornell Institute for High-Energy Phenomenology, Cornell University, Ithaca, NY 14853 

Abstract 

The introduction of T parity dramatically improves the consistency of Little Higgs 
models with precision electroweak data, and renders the lightest T-odd particle (LTP) 
stable. In the Littlest Higgs model with T parity, the LTP is typically the T-odd heavy 
photon, which is weakly interacting and can play the role of dark matter. We analyze 
the relic abundance of the heavy photon, including its coannihilations with other T- 
odd particles, and map out the regions of the parameter space where it can account for 
the observed dark matter. We evaluate the prospects for direct and indirect discovery 
of the heavy photon dark matter. The direct detection rates are quite low and a 
substantial improvement in experimental sensitivity would be required for observation. 
A substantial flux of energetic gamma rays is produced in the annihilation of the heavy 
photons in the galactic halo. This flux can be observed by the GLAST telescope, and, 
if the distribution of dark matter in the halo is favorable, by ground-based telescope 
arrays such as VERITAS and HESS. 



1 Introduction 



It has now been firmly established that about 25% of the energy density in the universe 
exists in the form of nonrelativistic, non-baryonic, non-luminous matter, so called "dark 
matter" [1]. The microscopic composition of dark matter remains a mystery, but it is clear 
that it cannot consist of any elementary particles that have been directly observed in the 
laboratory so far^. Many theories which extend the standard model (SM) of electroweak 
interactions contain new particles with the right properties to play the role of dark matter; 
perhaps the best known example is the lightest neutralino of supersymmetric (SUSY) models. 

Recently, a new class of theories extending the SM at the TeV scale, "Little Higgs" 
(LH) models, has been proposed [2] (for reviews, sec [3, 4]). The LH models contain a light 
(possibly composite) Higgs boson, as well as additional gauge bosons, fermions, and scalar 
particles at the TeV scale. The Higgs is a pscudo-Nambu-Goldstone boson, corresponding 
to a global symmetry spontaneously broken at a scale / ~ 1 TeV. The global symmetry 
is also broken explicitly by the gauge and Yukawa couphngs of the Higgs. As a result of 
this breaking, the Higgs acquires a potential; however, the leading (one-loop, quadratically 
divergent) contribution to this potential vanishes due to the special "collective" nature of the 
explicit global symmetry breaking, and the lightness of the Higgs can be achieved without 
fine-tuning. The dynamics of the Higgs and other degrees of freedom relevant at the TeV scale 
is described by a non-hnear sigma model (nlsm), vahd up to the cutoff scale A ~ Anf ~ 10 
TeV. In particular, the Higgs mass term is dominated by a one-loop, logarithmically enhanced 
contribution from the top sector, which can be computed within the nlsm and shown to have 
the correct sign to trigger electroweak symmetry breaking, providing a simple and attractive 
explanation of this phenomenon. Above the cutoff scale, the model needs to be embedded in 
a more fundamental theory; however, for many phenomenological applications, including the 
analysis of this paper, the details of that theory are not relevant and the nlsm description 
suffices. 

The Littlest Higgs model [2] is simple and economical, and it has been the focus of 
most phenomenological analyses to date [5]. Unfortunately, the model suffers from severe 
constraints from precision electroweak fits, due to the large corrections to low-energy ob- 
servables from the tree-level exchanges of the non-SM TeV-scale gauge bosons and the small 
but non- vanishing weak-triplet Higgs vacuum expectation value (vev) [6]. To alleviate this 
difficulty, the symmetry of the theory can be enhanced to include a Z2 discrete symmetry, 
named "T parity" [7]. In the Littlest Higgs model with T parity (LHT) [8], the non-SM 
gauge bosons and the triplet Higgs are T-odd, forbidding all tree-level corrections to preci- 
sion electroweak observables^. Loop corrections to precision electroweak observables in the 
LHT model were considered in [9], and the model was shown to give acceptable electroweak 

"'^In principle, it remains possible that dark matter consists of microscopic black holes made out of ordinary 
particles. However, we do not know of a compelling cosmological scenario in which this possibility is realized. 

^In the version of the model considered here, there is one non-SM T-even state, the "heavy top" T+. 
However it only contributes at tree level to observables involving the weak interactions of the top quark, 
which are at present unconstrained. 
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fits in large regions of parameter space compatible with naturalness. 

An interesting side effect of T parity is that the lightest T-odd particle (LTP) is guar- 
anteed to be stable. Analyzing the spectrum of the model, Hubisz and Meade [13] have 
argued that the LTP is likely to be the electricaly neutral, weakly interacting "heavy pho- 
ton" (or, more precisely, the T-odd partner of the hypercharge gauge boson) Bh- This 
particle is an attractive dark matter candidate, and initial calculations [13] showed that its 
relic abundance is within the observed range for reasonable choices of model parameters.^ 
In this paper, we will present a somewhat more detailed relic density calculation, including 
the possibility of coannihilations between the Bh and other T-odd particles. We will then 
discuss the prospects for direct and indirect detection of the heavy photon dark matter. 

2 The Model 

Our analysis will be performed within the framework of the Littlest Higgs model with T 
parity, which has recently been studied in Refs. [9, 13]. Let us briefly sketch the salient 
features of the model relevant here; for more details, see [9, 13] or the review article [4]. 

The model is based on an SU{5)/ S0{5) global symmetry breaking pattern; the Higgs 
doublet of the SM is identified with a subset of the Goldstone boson fields associated with 
this breaking. The symmetry breaking occurs at a scale / ~ 1 TeV. An [5*^7(2) x ^7(1)]^ 
subgroup of the SU (5) is gauged; this is broken at the scale / down to the diagonal subgroup, 
SU{2)l X f/(l)y, identified with the SM electroweak gauge group. The extended gauge 
structure results in four additional gauge bosons at the TeV scale, W^, Wfj and Bh-^ 

T parity is an automorphism which exchanges the [SU{2) x U{1)]^ and [SU{2) x C/(l)]2 
gauge fields; under this transformation, the TeV-scale gauge bosons are odd, whereas the 
SM gauge bosons are even. The odd gauge bosons have masses 

M{W%) ^ gf, M = M{Bh) ~ ^ ~ 0.16/, (1) 

where g and g' are the SM SU{2)l and U{1)y gauge couplings, and the normalization of 
/ is the same as in Ref. [9]. (Electroweak symmetry breaking at the scale v <^ f induces 
corrections to these formulas of order v"^ / p.) The "heavy photon" Bh is the lightest new 
gauge boson, and in fact is quite light compared to /. Since the masses of the other T-odd 
particles are generically of order /, we will assume that the Bh is the lightest T-odd particle 
(LTP), and it will play the role of dark matter candidate. The only direct coupling of the 
heavy photon to the SM sector is via the Higgs, resulting in weak-strength cross sections 
for Bh scattering into SM states. The heavy photon then provides yet another exphcit 

^While the LHT dark matter candidate is a spin-1 heavy photon, this is not an unambiguous prediction 

of Little Higgs models. For example, the "Simplest Little Higgs" models [10] supplemented by T-parity may 
contain a stable heavy neutrino which can play the role of dark matter [11], while closely related "theory 
space" models can give rise to a scalar WIMP dark matter candidate [12]. 

^The Wfj and Bh fields mix to form the two neutral mass eigenstates; however, the mixing angle is of 
order v/ f and can typically be neglected. 
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example of a weakly interacting massive particle (WIMP) dark matter candidate, and it is 
not surprising that wc will find reasonable regions of parameter space where it can account 
for all of the observed dark matter. For later convenience, we denote the mass of this particle 
by M. The range of the allowed values for this parameter is determined by the precision 
electroweak constraints, which put a lower bound on /, typically of about 600 GeV [9]. 
While there is no firm upper bound on /, we will assume / ^ 2 TeV to avoid reintroducing 
fine tuning in the Higgs sector. Using Eq. (1), this corresponds to the WIMP masses in the 
range 

100 GeV < M < 300 GeV. (2) 

In the scalar sector, the model contains an additional T-odd weak-triplet field 0, which has 
a mass of order / and no vacuum expectation value. In the fermion sector, each SM doublet 
{Q1 and Lj, where a = 1 ... 3 is a color index and i = 1 ... 3 is a generation index), acquires 
a T-odd partner, Q1 and Lj. The masses of these particles are also free parameters,^ with 
the natural scale set by /. To avoid proliferation of parameters, we will assume a universal 
T-odd fermion mass M for both lepton and quark partners; we will require M > M to avoid 
charged or colored LTPs, and assume M ^ 300 GeV, since otherwise the colored T-odd 
particles would have been detected in the squark searches at the Tevatron. In addition, 
non-observation of four-fermion operator corrections to SM processes such as e^e~ — > qq 
places an upper bound on the T-odd fermion masses [9] : 

MTeV<4.8/^,v, (3) 

where M and / are expressed in units of TeV. To cancel the one-loop quadratic divergence 
in the Higgs mass due to top loops, two additional new fermions are required in the top 
sector, the T-even T+ and the T-odd T_.® Their masses are related by 

so that there is just one additional independent parameter in this sector. We will choose 
it to be Mt_ , and assume Mt_ > M to avoid a charged LTP. The couplings of the heavy 
photon Bh which will be used in the calculations of this paper are summarized in Table 1. 



3 Relic Density Calculation 

In the early universe, the heavy photons are in equilibrium with the rest of the cosmic fiuid. In 
the simplest case of generic (non-degenerate) T-odd particle mass spectrum, the equilibrium 

^If the flavor structure of the T-odd quark mass matrix is generic, with order-one flavor mixing angles, 
the masses of the T-odd quarks need to be degenerate at the few per cent level [14] . 

^In Ref. [15], a variation of the model has been constructed where a single T-odd top partner is sufficient 
to cancel the divergences. Since the top sector will only play a minor role in the analysis of this paper, we 
expect our results to hold, at least qualitatively, in that model. 
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i (1) g'Y sin a (sin a^P^ + Pr) 



Table 1: Interaction vertices involving the heavy photon Bh that appear in the calculations 
of this paper. Here a = cos~ ^ {Mt_ / Mt^) , and Y = 1/10. 
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Figure 1: The leading 2^2 processes which maintain the heavy photon in equilibrium with 
the rest of the cosmic fluid at high temperatures. 



is maintained via the heavy photon pair-annihilation and pair-creation reactions; the leading 
2^2 processes that contribute are shown in Fig. 1. The present relic abundance of heavy 
photons is determined by the behavior of pair-annihilation rates in the non-relativistic limit, 
namely, by the sum of the quantities 



a(X) = limaiBnBH X)u, 



(5) 



over all possible final states X. Here, u is the relative velocity of the annihilating particles. 
Note that, unlike the bino-like neutralinos typically predicted by the constrained minimal 
supersymmetric standard model (cMSSM), the s-wave annihilation of the heavy photons is 
unsuppressed: in the language of Ref. [16], the heavy photons are "s annihilators" , analogous 
to the Kaluza-Klein photons of the "universal extra dimensions" (UED) model [17, 18]. It 
is straightforward to compute a{X) using the Feynman rules in Table 1. We obtain 



a{W^W^) 
a(ZZ) 



27ra2 




M' 




3 cos"^ 9w 


{AM' 


-ml)' 


+ mlTl 
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where /ij = mf/M'^, 9w is the SM weak mixing angle, and and Vh are the mass and the 
width of the SM Higgs boson. If M > rrit, the LTPs can also annihilate into pairs of top 
quarks; ignoring the contribution from the t- and it-channel T_ exchanges, we obtain^ 

If M > m/i, annihilation into a pair of Higgs bosons is possible, with the cross section 

fih{l + fJ'h/S) 1 



^ ' 2 COS^ Qyv 



(4M2 - m|)2 + m^r^ 24M4 



a/i - y^K. (8) 



Finally, LTPs can also annihilate into light SM fermions via t-channel exchanges of the T- 
odd fermions; this channel was not included in the analysis of Ref. [13]. For a fermion / 
(/ = z/, M, d) we obtain 

where — 1 for leptons and 3 for quarks, and Y — 1/10 is the Bffff couphng in units of g'. 
Because of the small value of Y, the annihilation into light fermions is strongly suppressed, 
even for relatively small values of M. The WMAP collaboration data [19] provides a precise 
determination of the present dark matter abundance: at two-sigma level, 

Qdmh^ = 0.111 ± 0.018. (10) 

For s annihilators, this translates into a determination of the quantity a = Y^x^^i^)'- — 
0.8 ± 0.1 pb. (The precise central value of a depends on the WIMP mass; however, this 
dependence is very mild, see Fig. 1 of Ref. [16].) Using this constraint and the above 
formulas, it is straightforward to map out the regions of the model parameter space where 
the heavy photons can account for all of the observed dark matter. The results are consistent 
with the updated analysis of Hubisz and Meade, see Fig. 3 of Ref. [13]. For given rrih, there 
are two values of M which result in the correct relic density. There is one solution on either 
side of the Higgs resonance. For WIMP masses in the interesting range, Eq. (2), these can 
be approximated by simple analytic expressions: 

mft 24 + 2.38M, or fa -83 + 1.89M, (11) 

where M and mh are in units of GeV. We will refer to these solutions as "low" and "high" , 
respectively. The analytic expressions (11) reproduce the values of M and nih consistent 
with the WMAP central value of ildm^^ with an error of at most a few GeV throughout the 
interesting parameter range. This accuracy will be sufficient for the analysis of detection 
prospects in Sections 4 and 5. 

^The r_ exchanges are negligible throughout most of the parameter space, but will nevertheless be fully 
included in the numerical calculation of the relic abundance described below. 
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Throughout the parameter space consistent with the WMAP value of the present dark 
matter density, the dominant heavy photon annihilation channels are W'^W~ and ZZ; the 
ti channel contributes at most about 5% of the total annihilation cross section, while the 
hh final state is always kinematically forbidden. Moreover, the ratio of the W^W~ and ZZ 
contributions is approximately 2:1, as is evident from Eqs. (6), so that a{W'^W~) 0.53 pb, 
a{ZZ) 0.27 pb throughout the parameter space. Since Bh is an s-annihilator, the same 
cross sections govern the rate of heavy photon annihilation in the galactic halo, which in 
turn determines the fluxes relevant for indirect detection, see Section 5. 

If some of the T-odd particles are approximately degenerate in mass with the heavy 
photon, the simple analysis above is no longer applicable, since coannihilation reactions 
between Bh and other states significantly affect the relic abundance. In the LHT model, 
the masses of the T-odd weak gauge bosons Wh and the triplet scalar are predicted 
unambiguously once the scale / and the Higgs mass ruh are fixed; these particles are always 
much heavier than the Bh and their effect is negligible. On the other hand, the common 
mass scale of the T-odd leptons and quarks M is a free parameter, and for M ~ M the 
coannihilations between these states and the Bh can be important. We have performed a 
more detailed analysis of the B^ relic density, taking this possibility into account. 

In the presence of coannihilations, the abundance calculation requires solving a system of 
coupled Boltzmann equations. We approached this problem numerically. The interactions of 
the LHT model were incorporated in the CalcHEP package [20], which was used to compute 
the scattering matrix elements for the appropriate processes. The rest of the calculation was 
performed using the DM++ package,*^ recently recently developed by one of us (AB). The pack- 
age first uses the matrix elements to compute the thermal averages (au), which determine 
the reaction rates entering the Boltzmann equation. Then, the freeze-out temperature of the 
dark matter is determined iteratively, using the Turner- Scherrer approximation [23]. Finally, 
the integral of (au) from freeze-out to present day (usually called J (xp) in the literature) 
is evaluated, providing the relic abundance. 

The results of this analysis are illustrated by Figure 2, which shows the contours of 
constant heavy photon relic density in the f — M (or, equivalently, M — M) plane. The 
typical situation for a heavy Higgs is shown in the left panel {rrih ~ 300 GeV). There are 
two regions in which the heavy photon can account for the observed dark matter: 

• The two vertical pair- annihilation bands, where the coannihilation processes are unim- 
portant. The heavy photon abundance in these regions is independent of M. The 
bands appear on either side of the s-channel Higgs resonance dominating the pair- 
annihilation processes, corresponding to the "high" and "low" solutions of Eq. (11). 
(The bands are analogous to the "Higgs funnel" region in the cMSSM.) 

^The DM++ package is currently being prepared for public release. It can be applied to compute the relic 
abundance of the WIMP for any particle physics model that can be incorporated in CalcHEP. The DM++ 
is inspired by the micrOMEGAs code [21], which was originally designed to compute the relic abundance of 
neutralinos in the MSSM. The recently developed new version of this code, micrOMEGAs2 . [22], is also 
applicable to any CalcHEP model defined by the user. This package is also being prepared for public release. 
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Figure 2: The contours of constant present abundance of the heavy photon LTP, Qltp^^, 
in the M — M plane. The Higgs mass is taken to be 300 GeV (left panel) and 120 GeV 
(right panel). The red and green contours correspond to the upper and lower bounds from 
WMAP, Eq. (10), assuming that the LTP makes up all of dark matter. The yellow and blue 
lines correspond to the LTP contributing 50% and 70%, respectively, of the measured dark 
matter density. The shaded region corresponds to a charged and/or colored LTP. 



• The coannihilation tail, where the heavy photon abundance is predominantly set by 
coannihilation processes. Since the T-odd fermions are assumed to be degenerate, all 
of them participate in the coannihilation reactions. The location and shape of this 
feature are similar to the tau coannihilation tail in cMSSM. 

As the Higgs mass is decreased, the pair-annihilation bands appear for lower WIMP 
masses, and for light Higgs (115-150 GeV) the "low" band disappears, since the required 
values of / are already ruled out by data. The "high" band persists until the Higgs mass 
is close to the current experimental bound. To illustrate this, consider the right panel of 
Fig. 2, where rrih = 120 GeV. The band between the two red lines (90 ^ M ^ 100 GeV) is 
allowed. Note that the behavior of the relic density as a function of M within this band is 
non-trivial: The relic density first drops with increasing M due to the fact that the threshold 
for the reaction BhBh ZZ is passed. It then bottoms out at a value consistent with the 
measured Vt^^Ji^, and begins increasing as increasing M further takes the center-of-mass 
energy away from the Higgs resonance, suppressing annihilation. Clearly, this situation is 
quite non-generic, and for somewhat higher irih the Z threshold becomes irrelevant and relic 
density is a uniformly increasing function of M in the "high" band. The coannihilation tail 
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Figure 3: The leading processes which contribute to the heavy photon-nucleon elastic scat- 
tering cross section relevant for direct dark matter detection experiments. 



is present for low as well as high values of m/j. The tail can be described by a simple analytic 
formula 

M^M + 20GeV, (12) 

which is approximately independent of the Higgs mass. 

It should be noted that the remaining free parameter of our model, the mass of the second 
T-odd top quark Mrp_, was fixed to be equal to /, so that Mt_ ^ M and this particle did 
not have an effect on the Bh relic abundance. We expect that a second coannihilation tail 
appears when Mt_ ~ M; the structure should be very similar to the one found above, with 
slight numerical differences due to smaller multiplicity of the coannihilating states. 



4 Direct Detection 

Direct dark matter detection experiments attempt to observe the recoil energy transfered 
to a target nucleus in an elastic collision with a WIMP. The null result of the current 
experiments places an upper bound on the cross section of elastic WIMP-nucleon scattering. 
In this section, we will discuss the implications of this bound for the LHT dark matter, and 
prospects for future discovery. 

The elastic scattering of the heavy photon on a nucleus receives contributions from several 
processes shown in Fig. 3. Consider first the scattering off gluons, which occurs via the Higgs 
exchange diagram (a). The Higgs-gluon coupling arises predominantly via a top quark loop, 
and has the form [24] 

^hgg = ■ ^^^^ 

where v = 246 GeV is the Higgs vev, and G^^, is the color field strength. The halo WIMPs 
are highly nonrelativistic (/3 ~ 10~^), and the momentum transfer in the reaction at hand 
is negligible compared to rrih. The WIMP-gluon interaction can then be described by an 
effective operator 

o; o; 1 

7 J^Ha^S'^LG'"''''. (14) 

6 COS'' 0]y ml ^ 
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In the chiral limit, the matrix element {n\G^\n) can be related to the nucleon mass nin [24], 
leading to an effective WIMP-nucleon vertex of the form 

2 

^eff = 07 %^ ^^Hai^g^n^n, (15) 

27 cos^ OwT^h 

where is the nucleon (neutron or proton) field. It is clear that this interaction only 
contributes to the spin-independent (SI) part of the WIMP-nucleon scattering cross section. 
Neglecting other contributions to the SI cross section (which, as we will argue below, are 
expected to be subdominant) , we obtain 

Ana'' mj 1 
729 cos^ Ow mi (M + rrinf ' ^ ' 

for both neutrons and protons. Since the scattering off nucleons in a given nucleus is coherent 
and the matrix elements for neutrons and protons are identical, the SI cross section for 
scattering off a nucleus of mass mjv is simply obtained from Eq. (16) by a substitution 

rrin ruN. 

The interaction of WIMPs with quarks is dominated by the T-odd quark exchange dia- 
grams, see Fig. 3 (b) and (c). (The Higgs exchange diagrams are suppressed due to small 
Yukawa couplings of quarks. In fact, it is well known that the Higgs-nucleon interaction 
is dominated by the Higgs-gluon coupling considered above.) The scattering amplitude is 
given by 



-i — tt;— <(P3)ei/(pi) u{p. 



C0S2 ew "'"^ ' [kl - M2 - M2 



Plu{p2), (17) 



where ki — pi + P2, ^2 = P2 — Ps- The qqEn coupling is flavor-independent, Y = 1/10, and 
the expression (17) is valid for every quark species. The amplitude contains two important 
physical scales: the weak scale, M ~ M ~ 100 GeV, and the QCD scale, Aqcd ^ 100 
MeV, which represents the typical energy and momentum of the quarks bound inside a 
stationary nucleus and, by a coincidence, the spatial momenta of the halo WIMPs: |p|i,3 ~ 
/3M ~ Aqcd- We will work to leading order in the ratio of these two scales. In this 
approximation, ki ^ —k2 ~ (M, 0), and the heavy photon polarization vectors are purely 
spatial, s^{pi^3) = (0,£i,3). The amplitude takes the form 

€ijke\eiu4-f''il - 7^)«2 , (18) 



cos2 9w M2 - M2 

corresponding to the coupling of the Bh spin with the vector and axial- vector quark currents. 
The axial current interaction corresponds to the coupling between the WIMP and quark 

spins, and gives rise to the spin- dependent (SD) part of the WIMP-nucleus scattering cross 
section. By the Wigner-Eckardt theorem, the quark axial current can be replaced by the 
nuclear spin operator s'^: 

'N\qYl'q\N) = 2s^^X,. (19) 
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For a nucleus of spin Jjv, the coefficients are given by 



A, = Ag,i^ + Ag„i|^, (20) 
Jn Jn 

where {Sp^n) /Jn is the fraction of the total nuclear spin carried by protons and neutrons, 
respectively, and the quantities Ag„ can be extracted from deep inelastic scattering data. 
We will use = Ad„ = 0.78 ± 0.02, = Am„ = -0.48 ± 0.02, As„ = Asp = 
—0.15 ± 0.02 [25]. The effective WIMP-nucleus spin-spin interaction can then be written as 

e,,fcSi,Si,*jv4*iv E K (21) 



yielding the SD cross section 



2 



- 3 cos^ (M + m^f (M2 - M^f ^^^^^ + 1 ^ ' " ^^^^ 



\q=u,d,s 



Now, consider the part of the amplitude (18) involving the quark vector current. Since 
the current is conserved, the contributions of each valence quark in a nucleon add coherently, 
and sea quarks do not contribute. The resulting WIMP-nucleon coupling is 

e,,fcS^Si,^„7'*n. (23) 



C0S2 ^iy(M2 -M2) 

This interaction is suppressed in the nonrelativistic limit, since ?7„7^-u„ ~ f^. In fact, it is of 
the same order as other contributions to the WlMP-quark scattering amplitude, suppressed 
by WIMP velocities or powers of Aqcd/^, which were neglected in our analysis. Therefore, 
its effect will be neglected. 

It should be noted that the SI interaction in Eq. (15), is parametrically suppressed with 
respect to the leading SD coupling, Eq. (21), by a factor of mj^/nih ~ Aqcd/^i and is 
formally of the same order as the contributions to the WIMP-quark interaction that were 
neglected in our analysis. Since the neglected terms contribute to the SI as well as SD 
interactions, one may question the validity of the SI cross section obtained in Eq. (16). 
Note, however, that the WIMP-quark interactions are additionally suppressed by a factor of 

= 0.01, not present in the WIMP-gluon couplings. Thus, while of the same order as (15) 
in terms of power counting, the neglected SI corrections from WIMP-quark interactions 
are expected to be numerically small. One interesting potential exception occurs in the 
coannihilation region, where the Y suppression could be compensated by the factor of M — 
M -C M in the propagator, and the WIMP-quark interactions could provide a significant 
correction to Eq. (16). A detailed analysis of this issue is reserved for future study. 

The SI elastic WIMP-nucleon scattering cross sections expected in the LHT models are 
plotted in Fig. 4, along with the current bound from the CDMS collaboration [26] (sohd red 
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Figure 4: The spin-independent (SI) WIMP-nucleon elastic scattering cross section in the 
pair-annihilation bands (left panel) and in the coannihilation region, for two values of rrih, 
120 and 300 GeV (right panel). The present [26] and projected [27] sensitivities of the CDMS 
experiment are also shown. 



lines) and the projected future sensitivity of SuperCDMS, stage C [27] (dashed red lines). 
We assume that the heavy photons account for all of the observed dark matter, and the two 
panels correspond to the two regions of parameter space which satisfy this constraint. The 
left panel shows the cross section expected in the pair-annihilation bands, with the two lines 
corresponding to the high and low solutions in Eq. (11). The right panel shows the cross 
section expected in the coannihilation tail for two values of the Higgs mass, 120 GeV and 
300 GeV. The two lines can be thought of as the upper and lower bounds on the expected 
cross section.^ While the predicted cross sections are two-three orders of magnitude below 
the present sensitivity, the expected improvements of the CDMS experiments will allow it to 
begin probing the interesting regions of the model parameter space in both pair-annihilation 
and coannihilation regions. 

Fig. 5 shows the spin-dependent cross sections predicted by the LHT model, along with 
the current bound from the NAIAD experiment [28] and its projected sensitivity [29]. In 
the pair-annihilation bands, the scale M is allowed to vary between 350 GeV and the upper 
bound given in Eq. (3). (Recall that for a given value of M, the scale / is fixed unambiguously 
by Eq. (1).) Unfortunately, the predicted SD cross sections are several orders of magnitude 
below the NAIAD sensitivity. 

^Note, however, that in the LHT model a heavy Higgs, ruh > 300 GeV, may be consistent with precision 
electroweak data in certain regions of parameter space where its contibution to the T parameter is partially 
cancelled by new physics contributions [9]. A heavier Higgs corresponds to smaller SI cross section. 
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Figure 5: The spin-dependent (SD) WIMP-proton elastic scattering cross section in the 
pair-annihilation bands (left panel) and in the coannihilation region (right panel). In the 
pair-annihilation bands, the scale M is allowed to vary between 350 GeV and the upper 
bound given in Eq. (3). The present [28] and projected [29] sensitivities of the NAIAD 
experiment are also shown. 



5 Indirect Detection via Anomalous Gamma Rays 

As discussed in Section 3, WIMP annihilation processes have to occur with approximately 
weak-scale cross sections to ensure that the relic abundance of WIMPs is consistent with 
observations. Since the heavy photons of the LHT model are s-annihilators, their annihilation 
rates are approximately velocity-independent in the nonrelativistic regime. This implies that 
the WIMPs collected, for example, in galactic halos, have a substantial probability to pair- 
annihilate, resulting in anomalous high-energy cosmic rays which could be distinguished from 
astrophysical backgrounds. In particular, high-energy gamma rays (photons) and positrons 
are considered to be the most promising experimental signatures. The gamma ray signal 
is particularly interesting because the gamma rays in the relevant energy range travel over 
galactic scales with no scattering, so that if the signal is observed, information about the 
WIMP (e.g. its mass) could be extracted from the spectrum. In this section, we will compute 
the gamma ray fluxes predicted by the LHT model, and evaluate their observability.^'' 

There are three principal mechanisms by which hard photons can be produced in WIMP 
annihilation: 

"'^^Positron fluxes from the heavy photon dark matter annihilation in the LHT model were recently con- 
sidered in Ref. [30]. 
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Figure 6: The diagrams which dominate the monochromatic photon pair-production in the 
Bh annihilation in the galactic halo. 



• Monochromatic photons produced via direct annihilation into a two body final state 
(77, /17 or Z7); 

• Photons radiated in the process of hadronization and fragmentation of strongly inter- 
acting particles produced either directly in WIMP annihilation (e.g. BhBh — > qq) or 
in hadronic decays of the primary annihilation products (e.g. BhBh —>■ ZZ followed 
by Z ^ qq); 

• Photons produced via radiation from a final state charged particle. 

Let us consider each of these mechanisms in turn in the LHT model. 

WIMPs being electrically neutral, production of monochromatic photons can only occur 
at loop level. In this paper, we will concentrate on the 77 final state, postponing the analysis 
of the Zj and hj channels for future work. (The photons produced in these reactions are 
separated in energy from the 77 photons due to non-zero masses of the Z and the Higgs.) 
The BhBh —>■ 77 process is dominated by the one-loop diagrams inducing the effective 
/177 vertex, see Fig. 6.^^ The corresponding cross section can be easily evaluated using the 
well-known formulas for the Higgs boson partial widths: 

a,,u ^ a (BhBh - 77) « = ^ {s - miy + 7^ ' ^^'^ 

where u is the relative velocity of the annihilating WIMPs, and s ~ 4M^ in the non- 
relativistic regime relevant for the galactic WIMP annihilation. The hat on P indicates that 
the substitution ruh ^/s should be performed in the standard expressions for on-shell 
Higgs decays [32, 33], and the loops of new particles present in the LHT model should be 
included. We obtain 

«^o^ s''/^ 2 

r'"-^'^'> = io2l^<;l-^- + -^'/^ + -^°l- 

^^A complete calculation would also include the contribution of the box diagrams with T-odd and T-even 
quarks running in the loop, analogous to the quark/squark boxes entering in the case of MSSM neutralino 
annihilation [31]. In the LHT case, this contribution is expected to be subdominant since the matrix element 
contains a factor of — 0.01. 



13 



where As denotes the contribution from loops of particles of spin s. These contributions are 
given by 



A 



1/2 



Ao 



E 



V2Qjyi 



.1/2 
W 



.1/2 



- 9 r. 



Fi/2{n); 



(26) 



where the sums run over all the charged particles of a given spin, and implicitly include 
summations over colors and other quantum numbers where necessary. The particles in the 
sums have masses mj and electric charges (in units of the electron charge) Qi, their trilinear 
couplings to the Higgs boson are given by AjV, 2/i/\/2, and Cig MwT]'^'^ , for particles of spin 
0, 1/2, and 1, respectively. (With these normalization choices, Ci — 1 for the SM W^, and 
yj's are the usual Yukawas for the SM fermions.) We have also defined Tj = 4mf/s. The 
functions Fs{t) are given by 



Fi(t) = 2 + 3r + 3r(2-r)/(r) 
Fi/2(r) = -2r(l + (l-T)/(T)), 
Fo(t) = t(1-t/(t)). 



(27) 



where 



fir) 




if T < 1. 



(28) 



Using these expressions, we find that the contributions of the T-odd states are subdominant 
compared to the SM loops. The contributions of the T-odd fermion loops and the T-even 
heavy top loop are suppressed because their couphng to the Higgs is of order j p . The 
contributions of charged T-odd heavy gauge bosons and scalars are suppressed due to their 
large masses, of order /. The deviation of the effective /i77 coupling from its Standard Model 
value due to these states is of order a few per cent.^^ Given the much larger astrophysical 
uncertanties inherent in the anomalous photon flux predictions, we will ignore these effects 
in our analysis. 

The monochromatic flux due to the 77 flnal state, observed by a telescope with a line of 
sight parametrized by ^ = (^, (^) and a fleld of view AO can be written as [35] 



1.1 X IQ-V 



-icm-2 



100 GeVy 



J(^, An)AO. 



(29) 



i^Thc deviations of the h ■ 
in detail in Ref. [34]. 



77 and gg 



^1 pby V M 
h vertices from the SM in the LHT model were recently analyzed 
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Figure 7: The flux of tlie monocliromatic pliotons from tlie reaction BhBh — > 77 in tfie 
pair-anniliilation bands (left panel) and the coannihilation region (right panel). On the left 
panel, the blue/upper and the orange/lower lines correspond to the high and low solutions 
in Eq. (11), respectively. On the right panel, the blue/upper line corresponds to rrih = 300 
GeV, the orange/lower line to rrih = 120 GeV. The plots assume J(\E', Ai7) Afi = 1; all fluxes 
scale linearly with this parameter. 



The function J contains the dependence of the flux on the halo dark matter density distri- 
bution: ^ 

J(*,Afi) = \ — I ^ ) / dn f p^dl. (30) 

^ ' ' 8.5 kpc \^0.3 GeV/cmV AnJi^n ^ ^ 

where I is the distance from the observer along the line of sight. Many models of the galactic 
halo predict a sharp peak in the dark matter density in the neighborhood of the galactic 
center, making the line of sight towards the center the preferred one for WIMP searches. 
However, the features of the predicted peak are highly model-dependent, resulting in a large 
uncertainty in the predicted J. For example, at AQ = 10~^ sr, characteristic of ground-based 
Atmospheric Cerenkov Telescopes (ACTs), typical values of J range from 10^ for the NFW 
profile [37] to about 10^ for the profile of Moore et.al. [38], and can be further enhanced by 
a factor of up to 10^ due to the effects of adiabatic compression [39]. 

The monochromatic photon fluxes (assuming JAQ = 1) predicted by the LHT model in 
the parameter regions where the heavy photon accounts for all of the observed dark matter, 
are shown in Fig. 7. The left panel corresponds to the pair-annihilation bands, and the 
right panel to the coannihilation region. Searches for gamma rays from WIMP annihilation 
have to be able to distinguish them from the astrophysical background. In the case of the 

""^^Note, however, that a powerful point-Uke source of ultra high energy gamma rays has been recently 
detected in the galactic center region [36]. The energy spectrum of this source, smooth and extending out 
to at least a few TeV, makes its interpretation in terms of WIMP annihilation unlikely. Detection of the 
potential gamma flux from WIMP annihilation in the same spatial region is clearly made more difficult by 
the presence of the source. 
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Figure 8: The fragmentation photon flux for M = 150, 200, 250 GeV (green, blue and red 
lines, respectively), in the pair-annihilation bands. The plot assumes Jl'^, AQ)AQ = 1; all 
fluxes scale linearly with this parameter. 



monochromatic photons, the signal is concentrated in a single bin (the energy uncertainty of 
the telescopes is about 10%, much larger than the intrinsic line width), and the background 
can be effectively measured in the neighbouring bins and subtracted. In the relevant energy 
range, the flux sensitivity for ground-based Atmospheric Cherenkov Telescopes (ACTs) such 
as VERITAS [40] and HESS [41] is estimated to be around (1 - 5) x lO'^^ cm~\ec~\ 
whereas the sensitivity of the upcoming space-based telescope GLAST is limited by statistics 
at 10"^'' cm~^sec~^, assuming that 10 events are required to claim discovery [42]. It is clear 
that the monochromatic flux predicted by the LHT model is beyond the reach of GLAST, 
but could be observed at the ACTs if the dark matter distribution in the halo exhibits a 
substantial spike or strong clumping, J ^ 10^ at AQ ^ 10~^. 

Let us now consider the component of the photon flux due to hadronization and frag- 
mentation of quarks produced in WIMP annihilation. As discussed in Section 3, the heavy 
photons predominantly annihilate into W and Z pairs; each of the vector bosons can in turn 
decay into a quark pair. The resulting photon spectra depend only on the initial energies 
of the W^s and Z's, and not on the details of the WIMP annihilation process. The spectra 
have been studied using PYTHIA (in the MSSM context), and a simple analytic fit has been 
presented in Ref. [35]: 

dN.y 0.73 _7 8^ 
ax x^-" 

where x = E^/M. This approximation is valid for both W~^W~ and ZZ final states. In the 
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100 GeV 
M 



) 



3 



J(*,AO)AO, (32) 



where we used the rehc density constraint, a{W'^W~) + a{ZZ) Ri 0.8 pb. The flux in the 
coannihilation region is much smaller. 

The fluxes predicted by Eq. (32) for several values of M are plotted in Fig. 8. The 
GLAST telescope is statistics-limited at energies above about 2 GeV, and would observe tens 
of events in this energy range for the heavy photon mass in the preferred range, assuming 
J(Af2)AQ = 1. (The flux scales linearly with this parameter combination.) One should 
keep in mind, however, that while the prospects for observing this signal are good, ruling 
out its interpretation in terms of conventional astrophysics could be challenging given the 
smooth, featureless nature of the fragmentation spectrum. Detailed studies of the angular 
distribution of these photons, in particular outside the galactic disk, will be needed. 

The ACTs have a higher energy threshold, typically about 50 GeV, and suffer from an 
irreducible background from electron-induced showers, about 10~^^ — 10~^° cm~^s~^GeV~^ 
in the relevant energy range (50 . . . 200 GeV) for Afi ~ 10"^ Using the extraploation 
of Ref. [35] to estimate the background, we find that the typical signal/background ratio 
expected at the ACTs, assuming Afl ~ 10""^ and J{Afl)Afl = 1, is only about 10"^. An 
observation of the fragmentation flux at the ACTs appears quite challenging, unless dark 
matter is strongly clustered at the galactic center or clumped. 

The third and final component of the gamma-ray flux from WIMP annihilation is the final 
state radiation (FSR) photons. The FSR fiux generally provides a robust signature of WIMP 
annihilation: it exists whenever the WIMPs have a sizable annihilation cross section into 
any charged states. The FSR photons have a continuous spectrum, in analogy to the quark 
fragmentation photons consdered above. In fact, at low energies, the fragmentation flux 
dominantes over the FSR component (unless WIMPs annihilate into purely leptonic states). 
At energies close to the WIMP mass, however, the fragmentation fiux drops sharply, and the 
FSR component typically dominates [43]. This is particularly interesting because the FSR 
spectrum typically possesses a sharp edge feature, abruptly dropping to zero at the maximal 
photon energy allowed by kinematics. The edge feature could help the experiments to discern 
this flux on top of the (a priori highly uncertain) astrophysical background, and provide 
a measurement of the WIMP mass [43]. In the LHT model, the dominant charged two- 
body annihilation channel is W~^W~, and correspondingly the reaction BhBh W~^W~'y 
provides the most important component of the FSR photon flux. The differential cross 
section for this process is given by 



BhBh W+W-'y) = a {BhB_ 



H 



W+W-) ^(x;/x^). 



(33) 



where x — lE^j ^ E^/M, //^ = (mw/My and 



a 



1 



1 



X (2x - 2 /x) log 



2(1 



x) — II — 2J {1 — x){l — X — iJ,) 



TT ^yl — jJi X 
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Figure 9: The FSR photon flux for M = 125, 150,200,250 GeV (left to right), in the pair- 
annihilation bands. (The fluxes for "high" and "low" solutions are essentially identical.) 
The plot assumes J(\E', Af2) Af2 = 1; all fluxes scale linearly with this parameter. 




1 -x)(l -x-/i)] , (34) 



V4-4/i + 3/i 

for < X < 1 — /i and for 1 — /i<a;<l. In the limit of large heavy photon mass, s ^ M, 
this expression reduces to 



-2 



2a; 1 — X 

TT X 



jr^x) = —^-1 log + - 1 + 0{fi)] . (35) 



The leading (logarithmically enhanced) term agrees with the result obtained in Ref. [43] 
using the Goldstone boson equivalence theorem. Note that the form of the photon emission 
factor JF, even in the large-s limit, depends on the theory being considered and on the 
initial state. For example, the XiXi ~^ W~^W~'j cross section in the MSSM, computed in 
Ref. [44], has a different leading logarithm behavior; this is related to the fact that the W 
bosons effectively become massless in this limit, inducing new infrared singularities. Thus, 
even though we chose to write the cross section (33) in a "factorized" form, there is no true 
factorization in the WW'j final state, in contrast to the ffj final states [43]. 
The flux of the FSR photons is given by 

§ = (5.6 X lo-v'cm-Gev-) (=^2^) {^^)\m,mm 

^ (36) 

where a(W^"*"iy ) is given in Eq. (6). Flux predictions for several representative values of M 
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are shown in Fig. 9. In each case, the flux drops abruptly at the maximal photon energy, 



El^ = (37) 

If this edge feature is observed, it would provide a robust signature of heavy photon annihila- 
tion, as well as a measurement of its mass. Note that the FSR and fragmentation components 
of the flux are comparable near the edge, so that the fractional drop in the "signal" flux at 
the edge is significant. Just as for the monochromatic and fragmentation photons, the sensi- 
tivity to the FSR flux at ACTs is limited by the background from electron-induced showers. 
Assuming a 10% uncertainty on the total flux measurement, the drop in the total flux as- 
sociated with the edge feature of the FSR spectrum can be observed if J ^ 10^ — 10^. The 
sensitivity at GLAST is limited by statistics, and an observation of the FSR edge requires 
even higher values of J. 

To summarize, we found that the best prospects for a discovery of anomalous gamma 
rays due to heavy photon annihilation in the Milky Way are offered by the GLAST telescope, 
which should be able to observe tens of fragmentation photons in the multi-GeV energy range. 
The fluxes of monochromatic and FSR photons, whose spectra would provide clear signatures 
for galactic WIMP annihilation (a bump and an edge, respectively), are significantly smaller. 
The prospects for their detection depend on the assumed halo profile; an observation by the 
ACTs such as VERITAS and HESS is possible if the dark matter density has a sharp peak 
at the galactic center or is strongly clumped, J ^ 10^ — 10^. 



6 Conclusions 

Little Higgs models provide an interesting alternative scenario for physics at the TeV scale, 
with a simple and attractive mechanism of radiative electroweak symmetry breaking. Many 
realistic models implementing the Little Higgs mechanism have been proposed; however, 
generically these models are ruled out by precision electroweak data, unless the scale / is 
in a few- TeV range which reintroduces fine-tuning. Little Higgs models with T parity avoid 
this difficulty. In this paper, we focused on the Littlest Higgs model with T parity (LHT), 
one of the simplest models in this class. T parity makes the lightest of the T-odd particles, 
the LTP, stable, enabling it to have a substantial abundance in today's universe in spite of its 
weak-scale mass. In the LHT model, the LTP is typically the heavy photon Bh, which can 
play the role of WIMP dark matter. We have computed the relic abundance of this particle, 
including coannihilation effects, and mapped out the regions of the parameter space where it 
has the correct relic abundance to account for all, or a substantial part, of the observed dark 
matter. These regions can be divided into the pair- annihilation bands, where the abundance 
is set by the Bh pair annihilation via s-channel Higgs resonance, and the coannihilation tail, 
where coannihilations of Bh with T-odd quarks Q and leptons L play the dominant role. 

In the second part of the paper, we evaluated the prospects for observing the heavy 
photon dark matter of the LHT model using direct and indirect detection techniques. Direct 
detection is quite difficult, due to the fact that the heavy photon predominantly couples to 
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the Standard Model states via the Higgs boson whose interactions with nucleons are weak. 
The elastic cross section of the Bh scattering on a nucleus in the region of parameter space 
consistent with the relic density constraint was found to be several orders of magnitude below 
the current sensitivity of direct detection searches such as CDMS. For indirect detection, we 
concentrated on the anomalous high-energy gamma ray signature. The predicted gamma ray 
flux depends sensitively on the distribution of dark matter in the halo. The best discovery 
prospect is offered by the GLAST telescope, which can observe the photons arising from 
the fragmentation of the W/Z bosons produced in the heavy photon annihilation. If dark 
matter distribution in the halo is favorable (in particular if it ehxibits a sharp spike near 
the galactic center, or is highly clumpy on short scales), ground-based telescopes such as 
VERITAS and HESS may also be able to observe a gamma ray signal. In this case, it might 
also be possible to observe the monochromatic and the FSR components of the photon flux, 
whose spectra exhibit well-defined features (a line and an edge, respectively) and would 
provide a smoking-gun evidence for the WIMP-related nature of the signal. 
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